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The scaling properties at low x of the proton DIS cross section and its charm component are 
analyzed with the help of the quality factor method. Scaling properties are tested both in the deep 
inelastic scattering data and in the structure functions reconstructed from CTEQ, MRST and GRV 
parametrisations of parton density functions. The results for DIS cross sections are fully compatible 
between data and parametrisations. Even with larger error bars, the charm component data favors 
the same geometric scaling properties as the ones of inclusive DIS. This is not the case for all 
parametrisations of the charm component. 

o 

o 
O 

Geometric scaling fli] is a remarkable empirical property verified by data on high energy deep inelastic scattering 
(DIS) i.e. virtual photon-proton cross-sections. One can indeed represent with reasonable accuracy the cross section 
a 1 p in the low Bjorken x regime by the formula 

<P'*(x,Q)=<t<'(t) , where r = log (gj^ry) (1) 

t is the scaling variable, Q the virtuality of the photon, and Y = log(l/x) the total rapidity in the 7*-proton system. 
I— ~~ '■ That empirical property is often considered as a hint for gluon saturation effects. Indeed, including non-linear effects 
due to high parton density in the evolution equations for parton distributions towards low x leads to such geometric 
scaling properties. In that context, Q S (Y) is called the saturation scale, and gives the typical scale for the onset 

' of nonlinear effects. Using the dipole factorization of DIS observables, one predicts geometric scaling to hold also 
for deeply virtual Compton scattering, exclusive vector meson production and inclusive diffraction, which has been 
verified empirically [2[ . Several types of geometric scaling or of generalizations of geometric scaling have been proposed 
in the literature, which correponds to various sets of approximations in the theoretical description of parton evolution 

■ and gluon saturation. 

t-H ' The quality factor (QF) introduced in |3[ is an interesting tool in order to discuss scaling properties. It allows 
to compare the validity of different scaling laws on a given data set, without any assumption concerning the shape 
of the scaling function. It has been used in 0, 0| in order to study the above-mentioned scaling properties in the 
J> , data for DIS observables. A standard fit to data tests locally in each bin the value of an observable, whereas the 
QF only considers point-to-point correlations for the observable. Hence, the QF method could give complementary 
informations, difficult to catch with a standard fit if the data are not very precise. For that reason, it is interesting to 
test with the QF method to what extend global fits to parton distribution functions (PDF) leads to the same scaling 
properties as the data. Many efforts have been made in the last few years by the MRST and CTEQ groups in order 
to implement heavy quark mass effects in a consistent way in their global fits. That justifies the comparison of data 
and PDF sets not only for the inclusive DIS cross section a ~ F2/Q 2 but also for its charm component a c ~ F^/Q 2 . 
The data for the bottom component F\fQ 2 are not suitable for a scaling analysis, since they contain too few points. 

The plan of the paper is as follows. In section II we review recent improvements in the treatment of heavy flavors 
in PDF's global fits, and we present the PDF sets used in the present study. In section III we describe the various 
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scaling hypotheses to be tested and their theoretical motivation. In section IV we present the Quality Factor method 
we shall use in the analysis and the set of data we consider for the fitting procedure. We provide and comment our 
results for F2 in section V and for F$ in section VI. The section VII is devoted to the summary and the discussion 
of the results. 



Let us briefly discuss the recent progresses concerning the treatment of heavy quark mass effects in PDFs global fits 1 . 
In the collinear factorization framework, one factorizes observables {e.g. proton structure functions) as a convolution 
of coefficient functions, describing partonic matrix elements, and PDFs for active flavors. In that formalism, active 
flavors are massless. Thus it should be relevant only if the quark masses corresponding to active flavors are small 
compared to the factorization scale usually chosen to be a hard scale of the problem, such as Q for DIS. The 
most standard pQCD scheme for heavy flavors is the fixed flavor number scheme (FFNS), where heavy flavors are not 
considered as active, and are generated only by boson-gluon fusion. It gives reliable results for Q = fip of the order of 
the heavy quark mass rrihf. However, for Q 2 3> rn^j higher order terms in the coefficient functions are enhanced by 
powers of log(Q 2 /m^), so that perturbation theory breaks down. In order to overcome this difficulty, variable flavor 
number schemes (VFNS) have been introduced. They consist in a succession of FFNS with one additional active 
flavor each time Q becomes larger than the mass of a heavy quark flavor. Two successive FFNS are matched with 
each other at the value Q 2 ~ m hf chosen for the change of scheme, in order to ensure that renormalization is done 
in a consistent way for all FFNS. In the the first implementations of VFNS, called zero-mass-VFNS (ZM-VFNS), 
the switching from one FFNS to another at Q 2 ~ m hf' wn ich is an unphysical change of factorization scheme, was 
coinciding with the threshold of production of the heavy flavor in the final state. That is not considered to be correct, 
and the latter should correspond to a threshold in the total energy W rather than in Q. That difficulty and a few 
others are solved if one goes to the general-mass- VFNS (GM-VFNS), build in @. One of the ingredients used in the 
GM-VFNS is the replacement of x by the rescaled variable 



when a pair of heavy quarks HH is produced in the final state, in order to ensure the right kinematical threshold: 
for W as small as 2m/j/, \ goes to 1 and thus the cross section for heavy quark pair production vanishes. 

The recent global fits of PDFs of the CTEQ and MRST groups are using such GM-VFNS. However, the details 
of the implementation differs between the two groups pj, l8|. In particular, they use different conventions to name 
the order of the analysis (see e.g. Q), such as LO, NLO, etc. In our study, we use the following general purpose 
GM-VFNS PDFs sets: the NLO CTEQ6.6M [|, the NNLO MRST2006 [nj, and the NLO MRST2004 II]. We 
consider also the NNLO MRST2004 which only uses an approximate implementation of GM-VFNS, and the 
older GRV98 [lj] PDFs set based on a FFNS with 3 flavors. 

In all of the abovementioned parametrisations, the charm production in DIS appears only pertubatively and the 
charm and anticharm PDFs (if any) start from zero at an intial factorization scale below the charm mass. Alternatively, 
some theoretical models for the nucleon content feature a non-perturbative intrinsic charm component. The CTEQ 
group studied that possibility, and released the PDFs sets CTEQ6.6C1 to CTEQ6.6C4 [g[, which extend CTEQ6.6M by 
including various types of intrinsic charm component. The CTEQ6.6C1 and CTEQ6.6C2 parametrisations correspond 
respectively to a moderate and a strong intrinsic charm contribution of the form predicted by the BHPS model [l3T ]. 
which is non-negligible only at large x. By contrast, the CTEQ6.6C3 and CTEQ6.6C4 parametrisations correspond 
respectively to a moderate and a strong intrinsic charm contribution proportional to the light sea quark contribution 
at the initial scale for the global fit. They are not motivated by theoretical models. 



Let us sketch the theoretical motivation for the different forms of scaling (see TABLE I) proposed for deep-inelastic 
scattering at high energy. The evolution of unintegrated parton densities towards large rapidity Y at fixed k 2 -, is given 
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II. PDFS WITH HEAVY QUARK MASS EFFECTS 




(2) 



III. SCALING VARIABLES FROM SATURATION PHYSICS 



by the BFKL [14j equation, provided that these unintegrated parton densities are small. Generic solutions of that 
linear evolution can be written as a sum of elementary wave solutions, with the help of Mellin transform. Each of 
these elementary wave solutions have scaling properties, but generically not their superposition. However, the BFKL 
evolution in Y leads to larger and larger unintegrated gluon density due to soft gluon radiation. At some point, the 
partons are so packed that additional soft gluons are emitted collectively. Thus, the soft gluons radiation is reduced 
by destructive interferences and the evolution becomes nonlinear in domain typically given by kx < Q S {Y). The 
BFKL equation is then generalized e.g. by the BK [15| or JIMWLK [16| nonlinear equations. For a review about 
that gluon saturation phenomenon, see [17] . Such a nonlinear evolution equation has the following property. Even 
in the kinematical domain in which the nonlinear terms are negligible compared to the linear ones, the nonlinearity 
of the equation constrains the solution, acting as a dynamical boundary condition to the linear BFKL evolution 
fli ]. The interplay between that specific boundary condition and the BFKL kernel selects dynamically [l9[ a specific 
wave solution (the critical one) of the BFKL equation in a window k\ > Q S (Y), which determines the evolution 
of Q S (Y). In that kinematical range, the solution of the nonlinear evolution loses memory of the initial condition. 
Hence, one finds that the unintegrated gluon distribution scales with k^/Q s (Y) above the saturation scale. Due to 
the fey-factorization in the linear regime, that property result in a geometric scaling ([l} of the DIS cross-section, and 
of other observables. 

In a fixed coupling a s approximation, that mechanism is well understood, and leads to a saturation scale log Q S (Y) ~ 
XY. This corresponds to the geometric scaling variable FC in TABLE I. If one goes to the running coupling case, the 
mechanism leading to the scaling holds, but is analytically under control only at large enough Q 2 , due to asymptotic 
freedom. One can deduce only the large Y behavior of the saturation scale, which is logQ s (F) ~ X\/Y '. Thus, one 
gets the predicted scaling property at large Y and Q 2 (with logQ 2 oc y/Y). However, the extrapolation of that scaling 
property towards the finite Q 2 and Y phenomenologically relevant domain is not unique. We consider in particular two 
different scaling variables, which both give the theoretical asymptotic behavior. The first one is a genuine geometric 
scaling variable (as defined in fl}). It is called RCI in TABLE I. The second scaling variable [2(| compatible with the 
asymptotic prediction is given in TABLE I and called RCII. Its value for the parameter A should be the square of the 
one for the RCI scaling. In the recent years, the possible impact on saturation of fluctuations due e.g. to Pomeron 
loops has been extensively discussed. They lead to a random saturation scale, but event-by-event the geometric scaling 
property is preserved. In the fixed coupling approximation 2 , the distribution of log<3 2 (Y") is approximately gaussian, 
with a variance proportionnal to Y. Finally, one expect the cross section to scale with the diffusive scaling variable 
[ll| (DS in TABLE I) which is the original geometric scaling variable divided by the dispersion \[Y . 



IV. THE QUALITY FACTOR METHOD 

In this section, we remind briefly the definition of the Quality Factor (QF) which we use to compare quantitatively 
the different scaling variables (see Table [J). We already used this method to compare the scaling results for the 
proton structure function F 2 , the deeply virtual Compton scattering (DVCS), the diffractive structure function, and 
the exclusive vector meson production data measured at HERA [H . 

Given a set of data points (Q 2 ,x,a — a(Q 2 ,x)) and a parametric scaling variable r = t(Q 2 , Y = log(l/a;); A) we 
want to know whether the cross-section can be parametrised as a function of the variable r only. Since the function 
of r that describes the data is not known, the QF has to be defined independently of the form of that function. The 
QF allows to quantitatively describe whether pairs of a = <r(Q 2 ,x) and r = t(Q 2 ,Y = log(l/x); A) lie on a single 
curve for a given parameter A. 

For a set of points (ui, Vi), where itj's are ordered, we introduce QF as follows 0] 



QF(X) 



E 



(Vi - Vi-i) 2 



(ui ~ Ui-!) 2 + e 2 



(3) 



where e is a small constant that prevents the sum from being infinite in case of two points have the same value of u. 
According to this definition, the contribution to the sum in ® is large when two successive points are close in u and 
far in v. Therefore, a set of points lying close to a unique curve is expected to have larger QF (smaller sum in ([3])) 
compared to a situation where the points are more scattered. 



2 In the running coupling case, such a scenario is expected to hold in the high Y limit However, the running of the coupling seems 

to suppress the fluctuations effects in the accessible rapidity range at collider experiments [2Jj| . 





scaling 


r formula 


parameters 


FC 


"Fixed Coupling" 


log Q 2 ~ AY 


A 


RCI 


"Running Coupling I" 


logQ 2 - Av 7 ^ 


A 


RCII 


"Running Coupling II" 


log(q 2 /A 2 ) - A log( j; /Aa) 


A 

A =0.2 GeV 


DS 


"Diffusive Scaling" 


\og,{Q*/K->)-XY 


A 

A =1 GeV 



TABLE I: Scaling variables used in the fits to deep inelastic scattering data Q. 

Since the cross-section in data differs by orders of magnitude and r is more or less linear in log(Q 2 ) (see Table [TJ) , 
we decided to take Ui — Tj(A) and Vi — log(erj). This ensures that low Q 2 data points contribute to the QF with a 
similar weight as higher Q 2 data points. The set (u^Vi) has to be ordered in u before entering the QF formula. In 
order to stay independent of the exact values of u and v, all the values are rescaled so that < Ui,Vi < 1. All the 
QF's in this Letter are calculated with e = 0.01. 

In order to test a scaling law r and say whether the points (<7i,Tj) lie on a single line or not, we search for the 
parameter A that minimises the 1/QF variable. The minimum value of 1/QF is obtained using the MINUIT package. 
Given the maximum value of the QF (minimum of 1/QF), we are able to directly compare different scaling laws. 

In this Letter, we aim to extend the study performed in Ref. |4{ by considering the different Fi parametrisations 
and the proton charm structure function data F£ measured at HERA and given by different parametrisations. 

V. FITS TO INCLUSIVE DIS 

We first aim to test the scaling quality using parametrisations of the proton structure function F2 using the 
parametrisations from the CTEQ, MRSTW and GRV groups discussed in section ILT1 

In a previous paper [3], we tested the scaling properties of F2/Q 2 using the data available from the HI [2~j| . 
ZEUS [23|, NMC [26| and E665 [27| experiments. We follow the same kinematical cuts for the parametrisations as 
we used for the experimental data (4j. We only consider data with x < 10 -2 and Q 2 in the range [3; 150] GeV 2 . The 
points with x > 10~ 2 are removed from the data since valence quark densities dominate here, and the formalism of 
saturation cannot apply in this kinematical region. Similarly, the upper Q 2 cut is introduced, while the lower Q 2 cut 
ensures that we stay away from the soft QCD domain. In this kinematical domain, 217 data points are used. In order 
to be able to directly compare scaling properties in the CTEQ, MRST and GRV parametrisations to the experimental 
data, we choose the same 217 pairs of x and Q 2 as in data. 

The QF method is more likely to prefer the parametrisations over data, since there is no statistical fluctuation in 
the parametrisations, contrary to data. It becomes clear from the definition that the statistically scattered values 
of structure functions in the experimental data lead to worse QF than a smooth parametrisation. In this sense, we 
cannot use the value of the QF to compare the scaling properties between experimental data and parametrisations. 
However, the QF can be directly compared among different parametrisations and can also test the differences among 
scaling laws tested on one particular data set. 

Table [TT] shows the fit results for all data and the parametrisations. Similarly as in data, the MRST and GRV 
parametrisations give similar QF for fixed coupling, running coupling I and running couping II. On the contrary, CTEQ 
favours fixed coupling. All parametrisations disfavour diffusive scaling, and all A values are consistent among different 
data sets. We can also see that GRV leads to smaller values of QF than MRST and CTEQ, but all parametrisations 
lead to a good QF. Scaling curves and fit results for experimental data and the CTEQ6.6M parametrisation, as an 
example, are shown in figure [TJ 

We note that the CTEQ, MRST and GRV parametrisations of the data lead to good scalings, whereas they do not 
contain any saturation hypothesis. One may wonder if such result comes from the chosen form of the parametrisation 
at the initial scale Qo, or if DGLAP evolution itself leads to a scaling in the HERA kinematical range. 

VI. FITS TO THE CHARM COMPONENT OF DIS 

In a second part of the paper, we want to study the scaling properties of the charm component of the structure 
function F^/Q 2 , both in data and in parametrisations. 



A. Fits to F2/Q 2 in data and parametrisations 



First we test the scaling properties using experimental data. The requirements on the kinematical domain remain 
the same as in the previous section. Now the lower Q 2 > 3 GeV 2 cut also allows to remove a part of the charm mass 
effects. We use the charm F£ measurements from the HI, ZEUS and EMC experiments Only 25 data points lie 
in the kinematical region we want to analyse. 

Since the statistics in the data is low, the fit results are not precise. Nevertheless, they still lead to clear results that 
are comparable to F2 fits. Table 111 shows the results of the fits to the F% data for different scaling laws. The results 
are found similar between F2 and F£ (see Fig. 2, top left plot). All A parameters are similar except for Diffusive 
Scaling. As in the case of the F2 scaling analysis (4j, Fixed Coupling, Running Coupling I and Running Coupling II 
give similar values of QF, and Diffusive Scaling is disfavoured. 

It is now interesting to check whether the scaling properties are also observed using the F£ parametrisation of the 
data, as seen in the previous section for the proton structure functions F2 . The results are given in Table III and Fig. 
2. Fig. 2 displays the values of the QF for Fixed Coupling and different parametrisations. The value of A favoured by 
the GRV and MRST parametrisations has a tendency to be larger (about 1) than the value favoured in data (about 
0.3). It is worth noticing however that the data show a smaller peak in the QF distribution around 1 as well which 
is disfavoured, and the MRST parametrisation another peak towards 0.3. The CTEQ parametrisation is definitely 
closer to the results in data leading to a value of A close to 0.5 for the CTEQ6.6M NLO parametrisation, and to 0.4 
when one introduces an additional sea-like intrinsic charm contribution 3 (CTEQ6.6C3 and CTEQ6.6C4). The results 
for the different scalings are given in Table III, and the conclusion remains unchanged. 

Let us now summarize what we obtained comparing data and simulation of F£. First of all, the value of A favoured 
in data is compatible with the one obtained for F2, DVCS... Of course the value of A shows a large error bar since 
the number of data points is only 25, but it is quite striking that the same value as for F2 is found. Everything 
looks as if the charm quark was behaving like any other light quark from the point of view of scaling. New data 
on i*2 which will be published soon by the HI and ZEUS collaborations are of great interest. The MRST and GRV 
parametrisations favour different values of A while the CTEQ ones are compatible with data, especially when an 
instrinsic charm component is added. 

The next step of our study relies more in understanding in more detail the differences between the parametrisations, 
since we recall that they lead to a good description of F^ itself. 

B. Comparison between different F£ parametrisations 

Taking the same 25 points as in data is not sufficient to investigate differences among the parametrisations. In 
order to analyse the parametrisations in detail, we choose to use the same x and Q 2 as for F2 which leads to 217 
points. 

The results for the 217 points are given in Table IIVI for the CTEQ parametrisation. For the sake of clarity, we 
choose to perform the study for Fixed Coupling geometric scaling only — the differences among different scaling laws 
are not large and do not modify the conclusion of the study. We note that scaling is indeed obtained in the different 
CTEQ parametrisations, and the value of A which maximises QF is compatible with the value measured in data. 
The CTEQ6.6C4 parametrisation, which has a strong sea-like intrinsic charm contribution, gives the A value that is 
closest to A = 0.3 found in data and leads to a much better QF than the other parametrisations. It is also worth 
noticing that the values of A are quite stable for different Q 2 cuts. 

A new striking result appear with the MRST and GRV parametrisations when one uses 217 points. The MRST 
2006 NNLO and the MRST 2004 NLO and NNLO parametrisations, as well as the GRV98 one, do not show any 
scaling at all. No value of A that maximises QF is found within an acceptable range. We now study the behaviour 
of the MRST and GRV parametrisations after different cuts on Q 2 . Increasing the cut to 5 GeV 2 does not help, 
therefore we choose to study the effects of higher Q 2 cuts. Table fVl gives the fit results in three different Q 2 ranges, 
namely [10; 150], [15; 150] and [25; 150] GeV 2 . Scaling properties are not observed using the MRST parametrisations 
even in the Q 2 range [10; 150] GeV 2 . The values of QF are much smaller than the ones obtained with CTEQ (see 
Table EJ). However, geometric scaling appears at higher Q 2 and the value of A ~ 0.65. The situation for the GRV98 
parametrisation is found to be similar as for MRST, but leading to a higher value of A. 



3 CTEQ6.6C1 and CTEQ6.6C2 lead to the same results as for CTEQ6.6M since the intrinsic charm contribution only appears at high x 
and does not affect our study performed for x < 10 — 2 . 



data set 


FC 


RC I 


RC II 


DS 


data 


A=0.33 


A=1.84 


A=3.44 


A=0.36 




QF=1.63 


QF=1.62 


gF=1.69 


QF=1.44 


GRV98 


A=0.33 


A=1.87 


A=3.35 


A=0.37 




QF=5.5 


QF=5.2 


QF=4.8 


QF=3.7 


MRST 2004 NLO 


A=0.34 


A=1.72 


A=3.31 


A=0.31 




QF=11.3 


QF=9.6 


QF=7.2 


QF=4.4 


MRST 2004 NNLO 


A=0.34 


A=1.72 


A=3.30 


A=0.36 




QF=11.2 


gF=n.3 


QF=7.6 


QF=4.3 


CTEQ 6.6 M 


A=0.35 


A=1.84 


A=3.40 


A=0.38 




QF=12.4 


QF=6.6 


QF=6.4 


QF=2.7 


CTEQ 6.6 CI 


A=0.35 


A=1.84 


A=3.40 


A=0.38 




QF=12.3 


QF=6.5 


QF=6.4 


QF=2.7 


CTEQ 6.6 C2 


A=0.35 


A=1.84 


A=3.40 


A=0.37 




QF=12.\ 


QF=6.4 


QF=6.3 


QF=2.7 


CTEQ 6.6 C3 


A=0.33 


A=1.84 


A=3.33 


A=0.31 




QF=11.7 


QF=6.2 


QF=7.2 


QF=2.9 


CTEQ 6.6 C4 


A=0.33 


A=1.66 


A=3.11 


A=0.38 




qf=i3.i 


QF=9.0 


QF=6.9 


QF=2.4 



TABLE II: Values of QF and A parameters for F2 data and the different F2 parametrisations from CTEQ, GRV and MRST 
and for fixed coupling, running coupling I, running coupling II and diffusive scalings. We note that all scalings lead to a good 
QF except for diffusive scaling and that the values of the A parameters are found similar in data and in the parametrisations 

As a conclusion to this second study, it looks like the MRST, GRV and CTEQ parametrisations of F% show a 
different behaviour from the scaling point of view. While the CTEQ parametrisation shows a scaling even at lower 
Q 2 (which seems observed in data), the MRST parametrisation only shows scaling at higher Q 2 . Further more, the 
GRV parametrisation leads to a much higher value of A than the one favoured in data. This difference might be due 
to the different forms of the gluon distribution at the initial scale Q 2 , chosen by the CTEQ and MRST groups. It will 
be worth to test this further using new measurements of F% which should be available soon. 

Finally, let us mention that we have checked that substitution in the scaling formulae of x by the rescaled variable 
X (see Formula ^) used in the GM-VFNS has a negligible effect in the considered kinematical ranges and does not 
change our results and conclusions. 

VII. CONCLUSIONS 

We have investigated with the Quality Factor method the scaling properties of F2 / Q 2 and of its charm component 
F 2 /Q 2 , both in the data and in a selection of PDFs parametrisations. The parametrisations lead to the same scaling 
properties of F 2 /Q 2 as in data for all types of scaling variables suggested by gluon saturation theory. Thanks to the 
precision of the F2 data, the global fits can indeed catch precisely the evolution of F 2 . The F 2 HERA data contain 
few points and still show large error bars. The QF method shows that the Fg data favors precisely the same value of 
the parameter A for each type of scaling as the one found by the QF fit on F 2 /Q 2 and on DVCS data 0. The CTEQ 
parametrisation gives the closest scaling properties for F 2 /Q 2 to the data, and the MRST and GRV parametrisations 
lead to favoured values of A higher than in data. 

To study further the differences between the parametrisations, we used 217 points, at the same x and Q 2 values 
as for F 2 . The parametrisations of the MRST and GRV groups have no scaling behavior at all if one includes points 
with Q 2 in the range [3; 10] GeV 2 . Only the CTEQ (and especially the CTEQ6.6C4 parametrisation, which includes 
a strong sea-like intrinsic charm component) has a scaling behavior with a value of A close to data. Moreover, it is 
the only parametrisation which leads to a very good scaling of F^jQ 2 for Q 2 down to 3 GeV 2 . The fact that the 
F 2 /Q 2 and F 2 /Q 2 data shows the same scaling properties seems to indicate that the behavior of the charm quark 
distribution in the proton is closer to light sea quarks in data than in most parametrisations. 

Of course, our study rely to a large extent on the QF analysis of the F% data which are not very precise. Hence, new 



data set 


FC 


RC I 


RC II 


DS 


data 


A=0.34 
QF=1.2\ 


A=1.72 
QF=1.20 


A=3.18 
QF=l.n 


A=0.22 
QF=1.16 


GRV98 


A=1.05 
QF=28.4 


A=5.67 
QF=25.2 


A >10 


A=0.17 
QF=9.9 


MRST 2004 NLO 


A=0.96 
QF=I0.0 


A=5.24 
QF=7.8 


A >10 


A=0.16 
QF=11.3 


MRST 2006 NNLO 


A=0.98 
QF=8.3 


A=5.35 
QF=6.7 


A >10 


A=0.16 
QF=9.9 


CTEQ 6.6 M 


A=0.49 
QF=28.3 


A=2.58 
QF=33.3 


A=7.73 
QF=29.I 


A=0.56 
QF=36.9 


CTEQ 6.6 C3 


A=0.48 
QF=46.I 


A=2.42 
QF=48.7 


A=6.03 
QF=35.2 


A=0.51 
QF=46.5 


CTEQ 6.6 C4 


A=0.42 
QF=51.7 


A=2.00 
QF=54.9 


A=4.50 
QF=40.1 


A=0.31 
QF=48.2 



TABLE III: Values of QF and the A parameters using F% data measured at HERA and the different F§ parametrisations for 
fixed coupling, running coupling I, running coupling II and diffusive scalings. We note that the values of A are similar than the 
ones found with the proton structure function F2 and that diffusive scaling is disfavoured. 





CTEQ 6.6 M 


CTEQ 6.6 CI 


CTEQ 6.6 C2 


CTEQ 6.6 C3 


CTEQ 6.6 C4 


Q 2 > 3 


A=0.70 


A=0.70 


A=0.68 


A=0.56 


A=0.44 




QF=0.8 


QF=0.8 


QF=0.8 


QF=3.5 


QF=14.6 


Q 2 >5 


A=0.61 


A=0.61 


A=0.61 


A=0.52 


A=0.42 




QF=5.0 


QF=5.0 


QF=5.0 


QF=12.5 


QF=19.6 


Q 2 > 10 


A=0.56 


A=0.56 


A=0.56 


A=0.50 


A=0.43 




QF=f9.2 


QF=19.3 


QF=I9.4 


QF=22.2 


QF=21.2 



TABLE IV: Values of QF and of the A parameter for F 2 C data and the different CTEQ parametrisations for different Q 2 lower 
cuts for fixed coupling scaling. We note that the scaling is observed in the CTEQ parametrisations and the value of A is stable 
for different Q 2 cuts. 

F£ data with improved statistics would be welcome, in order to check accurately if the scaling properties of F2/Q 2 
and F2/Q 2 are indeed the same. However, as the coincidence for the optimal value of A for F£/Q 2 , F2/Q 2 and DVCS 
is very striking, let us consider that there could be a deep reason for that, and discuss possible explanations. 

Since the behavior of the GRV98, which is based on a FFNS, and of recent MRST parametrisations based on 
GM-VFNS is similar, we conclude that the type of scheme for heavy flavors does not change much our discussion. 
Similarly, the order (NLO or NNLO) of the global fit seems not to be essential for the scaling properties. Thus, the 
similarity of the scaling properties of F 2 /Q 2 and F^jQ 2 in the data seems to result from a theoretical ingredient at 
low Q 2 or low x, which is missing in all of the usual global fits. 

The parametrisation which gives the best scaling properties and is the closest to data is CTEQ6.6C4. It could be 
tempting to say that our study gives hints for a strong sea-like intrinsic charm component in the proton. However, 
this is not the only possible explanation. For example, low x resummations and/or saturation effects would lead 
to a larger and more stable gluon distribution at low x and low Q 2 in a global fit. This would enhance the charm 
contribution produced by boson gluon fusion in that kinematical range without the need of intrinsic charm. Moreover, 
in that scenario, the light and heavy sea quark distributions would be both driven by the gluon distribution rather 
than by their initial conditions, and thus would naturally behave in a similar way. 





MRST 2004 NLO 


MRST 2004 NNLO 


MRST 2006 NNLO 


GRV98 


Q 2 > 10 


A=0.80 


A=0.80 


A=0.80 


A=0.94 




QF=2.7 


QF=3.0 


QF=2.5 


QF=4.4 


Q 2 > 15 


A=0.68 


A=0.70 


A=0.71 


A=0.90 




QF=7.9 


QF=7.9 


QF=6.7 


QF=8.7 


Q 2 > 25 


A=0.66 


A=0.64 


A=0.67 


A=0.87 




QF=18.5 


QF=18.7 


QF=18.7 


QF=21.2 



TABLE V: Values of QF and of the A parameter for F$ and the different MRST and GRV parametrisations for different Q 2 
lower cuts for fixed coupling scaling. Scaling is not observed for a Q 2 cut lower than 10 GeV 2 in the parametrisations. We 
note that the quality of scaling at low Q 2 is poor. For a higher Q 2 cut, scaling is observed but the value of the A parameter is 
different from the one in data for the GRV parametrisation. 
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FIG. 1: Different scalings for F 2 /Q 2 measured data and the CTEQ6.6M parametrisation for Q 2 > 3 GeV 2 : 

Upper curves: Scaling curve for "Fixed Coupling" (left) and QF against A for different scalings laws in data (right). Lower 
curves: same studies for the CTEQ6.6M parametrisation. Results are found similar for the different CTEQ, GRV and MRST 
parametrisations and are similar to data. 
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FIG. 2: Quality factor as a function of A for Fixed Coupling using different data sets. Quality factor peaks at the same value 
in F2 and F 2 C data. Quality factor in the F 2 C parametrisations peaks at higher values of A than in F 2 C data. 



